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A  low-density  ammonium  nitrate/fuel  oil  (ANFO)  explosive  was  de¬ 
veloped  by  diluting  the  normal  density  explosive  with  low-density 
polystyrene  beads.  Stable  detonations  were  achieved  in  mixtures 
having  densities  over  the  range  0.5  to  0.9  gm/cm-*.  The  detonation 
velocity  was  found  to  decrease  approximately  linearly  with  de¬ 
creasing  density.  The  explosive  was  found  to  have  a  long  reaction 
time  and  required  large  charges  with  propagation  distances  on  the 
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tion.  The  buildup  of  the  detonation  wave  to  its  steady- 
state  velocity  was  successfully  modeled  using  a  grain¬ 
burning  reaction  rate  theory  with  a  reaction  time  constant 
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SECTION  1 
INTRODUCTION 


Data  from  the  MINE  THROW  I  Event  indicate  that  the  detona¬ 
tion  velocity  and  pressure  for  normal-density  ANFO  are  consider¬ 
ably  higher  than  was  expected  from  values  published  in  the 
literature.  Pressures  were  measured  on  the  order  of  90  to  100 
kbar,  which  is  consistent  with  predictions  made  using  an  empirical 
technique  (References  1  through  4)  discussed  later  in  this  report. 
This  has  important  consequences  for  the  design  of  the  MINE  THROW 
II  simulation  charge.  The  outer  surface  of  the  simulation  charge 
is  located  at  the  peak  stress  contour  calculated  .or  the  nuclear 
event  corresponding  to  the  detonation  pressure  of  the  explosive 
used  for  the  simulation.  Therefore,  if  normal-density  ANFO  were 
used,  it  would  be  necessary  to  reduce  the  radius  of  the  charge 
outer  contour  to  correspond  to  the  90  to  100  kbar  peak  stress  con¬ 
tour.  However,  the  volume  enclosed  by  this  contour  is  so  small 
that  sufficient  explosive  cannot  be  put  into  the  cavity  to  match 
the  total  impulse  along  that  contour.  The  situation  is  made  worse 
by  the  fact  that  at  this  smaller  radius  more  energy  is  coupled 
to  the  ground  by  the  nuclear  event,  necessitating  the  use  of  a 
larger  quantity  of  explosive  for  the  simulation. 

To  circumvent  this  problem  a  low-density  ammonium  nitrate/ 
fuel  oil  explosive  with  lower  detonation  velocity  and  pressure 
has  been  developed  for  use  in  the  MINE  1HROW  II  Event.  This  has 
been  achieved  by  diluting  standard  prilled  ANFO  with  low-density 
styrofoam  beads.  Because  of  thair  low  density,  the  beads 
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contribute  very  little  mass  (typically  less  than  1  percent)  to 

the  explosive,  and  therefore  should  not  significantly  affect  the 

detonation  reaction.  Stable  detonations  have  been  achieved  in 

3 

mixtures  with  a  density  as  low  as  0.5  gm/cm  (compared  with  a 
normal  ANFO  density  of  0.86  to  0.90  gm/cm3) . 

A  survey  was  made  of  materials  suitable  for  diluting  the 
ANFO,  and  three  candidate  materials  were  chosen.  These  are 
listed  below  with  their  approximate  density  and  cost. 


Diluent 

Styrofoam  beads 
Saran  microballoons 
Bakelite  microballoons 


Density 

1  lb/ft3 

2  lb/ft3 
20  lb/ft3 


Cost 

$  1.50/ft3 
$10. 00/ft3 
$30. 00/ft3 


The  three  candidate  diluents  were  tested  and  each  was  found  to 
have  serious  disadvantages.  The  Saran  microballoons  tend  to 
clump  together  and  therefore  are  difficult  to  mix  uniformly  with 
the  ANFO.  Also,  because  of  their  small  individual  size,  a  large 
volume  of  microballoons  must  be  added  to  significantly  reduce 
the  density  of  the  explosive.  A  large  volume  of  Bakelite  micro¬ 
balloons  was  also  required,  although  they  mixed  with  ANFO  more 
easily.  Because  of  the  relatively  high  density  of  the  Bakelite 
microballoons  a  non-negligible  mass  of  material  was  added  to  the 
explosive,  which  is  undesirable.  In  both  cases  the  costs  were 
prohibitively  high. 


Styrofoam  beads  have  several  desirable  features.  They  are 
low  in  density,  low  in  cost,  and  the  volume  and  mass  of  material 
necessary  to  lower  the  explosive  density  is  less  than  for  the 
other  two  diluents.  The  major  problem  is  that  the  styrofoam 
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beads  tend  to  segregate  as  the  mixture  is  agitated,  causing 
variations  in  density  throughout  the  mix.  However,  by  mixing 
tbs  correct  proportion  of  styrofoam  beads  into  each  individual 
beg  of  AMFO,  this  effect  can  be  minimized.  Although  the  density 
of  >  the  mixture  within  tl.3  bag  may  not  be  uniform,  the  variations 
rill  be  on  a  scale  that  is  small  compared  to  the  charge  size 
foi  MINE  THROW  II  and  should  not  produce  a  significant  effect. 
This  has  been  verified  in  the  test  program  described  in  this 
report. 
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SECTION  2 
TEST  PROGRAM 


2.1  SMALL-SCALE  FEASIBILITY  TESTS 

To  test  the  feasibility  of  using  a  diluent  to  produce  a 
low-density  ANFO  explosive,  a  series  of  small-scale  tests  was 
performed.  The  purpose  of  the  tests  was  to  determine  the 
detonobility  of  diluted  ANFO  mixtures  at  various  densities  and 
to  provide  preliminary  data  to  help  in  designing  larger  tests. 

The  experimental  configuration  is  shewn  in  Figure  1.  The 
explosive  was  initiated  with  a  1-pound  pentolite  charge  and  '•.he 
detonation  trajectory  along  the  centerline  of  the  charge  was 
measured  with  ionization  switches.  Styrofoam  beads  were  used 
as  the  diluent  in  all  tests  except  one  (pq  -  0.71)  in  which 
Saran  microballoons  were  used.  Densities  ranging  from  0.5  gm/cm 
up  to  a  normal  density  of  0.88  gm/cm^  were  tested  and  no  problem 
in  achieving  a  detonation  was  found.  The  sensitivity  of  the 
explosive  increased  with  a  decrease  in  density,  and  a  systemati¬ 
cally  lower  final-detonation  velocity  was  observed  for  the  lower 
density  mixtures.  Figure  2  shows  values  of  the  final  detonation 
velocity  measured  at  the  end  of  the  charge  after  the  detonation 
had  propagated  approximately  80  cm.  These  are  plotted  as  a 
function  of  charge  density  and  are  compared  with  an  empirical 
prediction  of  the  steady-state  detonation  velocity  of  ANFO 
(discussed  later  in  this  report).  Figure  3  shows  a  comparison 
of  detonation  trajectories  in  three  different  density  mixtures. 
It  is  clear  from  this  figure  that  the  detonation  velocity  is 
still  increasing  and  that  a  larger  charge  size  is  necessary  to 
determine  the  steady-state  detonation  velocities. 

Preceding  page  blank 
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Figure 


Primacor 


1  Experimental  configuration  for  small-scale  detonation 
tests  of  low  density  AMFO. 


2.2  LARGE-SCALE  TESTS 


To  determine  the  steady-state  detonation  characteristics  of 
diluted  ANFO  as  a  function  of  density ,  a  series  of  large-scale 
tests  were  conducted  at  the  Nevada  Test  Site.  The  test  configura¬ 
tion/  shown  in  Figure  4/  consisted  of  an  8-foot  cube  of  explosive 
detonated  on  one  side  by  an  array  of  pentolite  charges  simulta¬ 
neously  initiated  to  produce  a  planar  detonation  front.  The 
array  of  initiators  is  shown  in  Figure  5.  The  ANFO  charge  was 
constructed  from  bags  of  the  diluted  ANFO  mixture  of  a  size 
equivalent  to  a  50-pound  bag  of  undiluted  ANFO.  Figure  6  shows 
the  weight  of  ANFO  per  bag  as  a  function  of  average  charge  density 
of  the  five  tests  conducted.  Also  shown  is  the  volume  of  styro¬ 
foam  beads  added  per  bag  for  those  five  tests.  Figure  7  shows  a 
typical  bag  of  the  diluted  ANFO  mixture.  It  is  clear  that  a 
significant  fraction  of  the  styrofoam  beads  has  become  segregated 
at  the  surface  of  the  bag.  The  effect  of  this  incomplete  mixing 
can  be  observed  in  the  detonation  trajectories  and  will  be  dis¬ 
cussed  later.  The  detonation  velocity  was  measured  with  an  array 
of  ionization  switches  oriented  as  shown  in  Figure  4.  The 
switches  were  taped  to  a  sheet  of  Mylar  to  assure  proper  position¬ 
ing  (Figure  8)  and  the  entire  assembly  was  placed  in  the  charge 
as  it  was  being  stacked  (Figure  9) .  Signals  from  the  ionization 
switches  were  recorded  on  raster  oscilloscopes  and  from  those 
records  a  time  history  of  the  detonation  front  was  determined 
for  each  test.  Manganin  pressure  gauges  were  mounted  on  the 
surface  of  the  charges  opposite  the  initiators  in  an  attempt  to 
measure  the  detonation  pressure.  However/  on  all  five  tests  the 
gauges  broke  before  a  peak  stress  could  be  recorded,  and  no  data 
were  obtained.  A  new  gauge  design  must  be  developed  and  tested 
prior  to  MINE  THROW  II  to  insure  that  pressure  measurements  are 
obtained  for  that  test. 
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ANFO  bags 


Plywood  box 


Figure  4  Test  configuration  for  the  large-scale  detonation 
tests  of  diluted  ANFO. 


Weight  of  ANFO  per  bag,  lb 


Figure  7  Typical  bags  of  diluted  ANFO  mixture. 
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One  plausible  explanation  tor  gauge  failure  on  all  the  tests 
is  the  possible  existence  of  a  very  irregular,  nonplanar  detona¬ 
tion  weve  front  in  the  ANFO.  Because  of  the  heterogeneous 
nature  of  the  diluted  ANFO  irregular  wave  front  is  likely, 
and  snail  jets  may  even  form  at  the  front  as  a  result  of  the 
interaction  of  the  ANFO  pr  11s  and  the  polystyrene  beads.  As 
this  irregular  front  impinges  on  the  manganin  gauge  element, 
large  deformations  could  be  produced,  causing  failure. 

If  this  is  the  cause  of  the  gauge  failure,  it  can  be 
circumvented  by  placing  a  thicker  buffer  between  the  explosive 
and  the  gauge  element.  A  layer  of  plexiglass,  or  some  other 
suitable  material  on  the  order  of  1  to  2  inches  thick  covering 
the  front  surface  of  the  gauge,  will  allow  the  irregularities 
in  the  shock  front  to  smooth  out  before  reaching  the  gauge. 

The  detonation  trajectories  measured  in  these  tests  are 
shown  in  Figures  10  through  14.  In  all  cases,  the  detonation 
velocity  increases  with  propagation  distance  until  a  constant 
value  is  achieved  at  a  distance  somewhere  between  100  and  150  cm. 
The  trajectories  are  more  irregular  than  those  measured  in  the 
small-scale  tests;  the  irregularity  increases  with  increased 
dilution  of  the  explosive  and  exhibits  a  periodic  structure  that 
has  a  characteristic  length  about  the  same  as  the  length  of  a 
bag  of  explosive.  The  irregularities  in  the  detonation  trajec¬ 
tories  are  probably  caused  by  the  incomplete  mixing  of  the  ANFO 
and  the  styrofoam  beads  within  each  bag.  Although  the  effects 
are  probably  not  important,  to  prevent  any  future  problems  a 
better  mixing  technique  should  be  developed  for  the  MINE  THROW  II 
explosive. 
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Figure  10  Detonation  trajectory  for  undiluted  94/6  ANFO. 
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Figure  11  Detonation  trajectory  for  diluted  ANFO. 
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Detonation  trajectory  for  diluted  ANFO 
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Figure  15  shows  a  composite  graph  of  the  detonation  trajec¬ 
tories  for  the  five  large-scale  tests.  The  comparison  shows  a 
consistent  trend  toward  lower  detonation  velocities  throughout 
the  trajectory  for  lower  density  explosives.  From  these  data, 
a  mixture  having  a  density  on  the  order  of  0.7  to  0.75  gm/cm3 
was  selected  for  the  HIKE  THROW  II  Event. 

To  design  the  MINE  THROW  II  charge,  a  mathematical  model 
describing  the  detonation  behavior  of  the  low-density  ANFO  was 
developed  from  these  data  to  be  used  in  the  finite-difference 
calculations  of  the  charge  performance.  Details  of  the  develop¬ 
ment  of  this  model  and  the  underlying  theory  are  presented  below. 
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Figure  15  Detonation  trajectories  for  diluted  ANFO 


SECTION  3 

MODELING  OF  DETONATION  PROPERTIES 


3.1  CHAPMAN- JOUGUET  DETONATIONS 

A  detonation  is  generally  represented  as  a  strong  shock, 
which  suddenly  compresses  and  heats  the  unreacted  explosive, 
followed  by  a  reaction  region  within  which  the  explosive  energy 
is  released.  This  is  shown  schematically  in  Figure  16.  The 
reaction  takes  place  over  a  finite  time,  t,  within  a  reaction 
zone  of  width  a.  The  energy  released  within  this  zone  provides 
the  potential  to  maintain  the  shock  front. 

In  a  manner  similar  to  deriving  the  Hugoniot  relations  for 
a  steady  shock  wave  in  an  inert  material,  the  equations  express¬ 
ing  conservation  of  mass,  momentum,  and  energy  can  be  used  to 
relate  quantities  on  either  side  of  the  reaction  zone  for  a 
steady-state  detonation.  These  equations,  in  the  steady-state 
reference  frame  moving  with  the  detonation  front,  are: 


p(D 

-  U)  =  pqD 

(i) 

P  + 

p(D  -  u)2  =  PQ  +  PqD2 

(2) 

E  + 

Pv  +  i  (D  -  u)  2  =  Q  +  E  +  Pv  +  J  D2 

2  o  o  o  2 

(3) 

where  p  is  density,  D  is  the  detonation  velocity,  u  is  the 
material  velocity,  P  is  the  pressure,  v  is  specific  volume, 
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Figure  16 


Schematic  representation  of  detonation  wave 


E  is  specific  energy,  and  Q  is  the  chemical  potential  energy 
(heat  of  reaction) .  These  can  be  manipulated  to  give  the  follow* 
ing  equations  for  the  variables  behind  the  reaction  zone: 


E  =  E  +  i  (I*  +  P  )  ( v  -  v)  +  Q  (7) 

0  2  O  O 


The  final  expression  is  the  Rankine-Hugoniot  equation  for  a 
detonation  wave.  From  these  equations  and  an  equation  of  state 
for  the  detonation  products,  P  *  f(V,E)  ,  the  Hugoniot.  curve  for 
the  detonation  wave  can  be  defined.  This  is  shown  in  Figure  17. 
It  can  be  shown  (Reference  5)  that  a  stable  steady-state  detona¬ 
tion  corresponds  to  a  unique  state  (P^,  v^)  on  the  detonation 
Hugoniot  for  which  the  detonation  wave  speed  D  is  equal  to  the 
sound  speed  (relative  to  the  laboratory  frame)  in  the  reaction 
products.  This  condition 


D  ■  u  +  C  (8) 

was  proposed  independently  by  Chapman  and  Jouguet  and  is  known 
as  *-he  Chapman-Jouguet  condition.  The  state  P^,  v^  is  known  as 
the  Chapman-Jouguet  state.  An  equivalent  condition  is  that  the 
Rayleigh  line  for  the  detonation  wave  be  tangent  to  the  Hugoniot 
of  the  products  as  shown  in  Figure  17  (Reference  5) . 


31 


Figure  17  Schematic  representation  of  a  detonation  in  the 
pressure- volume  plane. 


If  the  equation  of  state  of  the  products  is  assumed  to  be 
that  of  a  polytropic  gas 


(y  -  1)  pE 


(9) 


then  Equations  4  through  9  can  be  used  to  define  the  variables 
at  the  Chapman- Jouguet  state  (ignoring  PQ  relative  to  P^) : 


2  <  Y  “  l)pQQ 


PqD 


Y  +  1 


(10) 


2(y  -  1)Q 


(11) 


u. 


Tv""- IT 

Y  +  l 


Q  = 


D 

~ T 


(12) 


■  i  -j'  v 
Y  +  1  o 


(13) 


Y  +  1 


(14) 


These  equations  are  not  restricted  to  a  polytropic  gas  descrip 
tion  of  the  products,  but  hold  in  general  when  y  is  defined  as 

/a  In  p\  „ 

YCJ  “  In  vys  (1 

determined  at  the  Chapman-Jouguet  state. 
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In  a  detonation  wave  the  Chapman-Jouguet  state  is  not 
reached  instantaneously  but  after  a  reaction  tine  x.  According 
to  an  idealized  theory,  the  unreacted  explosive  is  shocked 
initially  to  sene  state  P*,  v*  shown  in  Figure  17.  The  chemical 
reaction  takes  place  over  a  time  x  with  the  products  reaching 
state  P^,  v^  as  the  reaction  goes  to  completion.  The  state 
P*,v*  corresponds  to  the  Von  Neumann  spike  at  the  detonation 
front  which  has  never  been  observed  experimentally.  In  real 
detonations,  that  state  may  never  be  reached  because  of  dissipa¬ 
tive  effects  (e.g.,  viscosity  and  heat  conduction)  in  the  shock 
front  (Reference  6) . 

For  most  conventional  military  explosives  the  chemical 
reaction  takes  place  very  quickly.  Since  in  most  applications 
times  of  interest  are  generally  long  compared  to  the  reaction 
time  and  the  charge  sizes  are  large  compared  to  the  reaction 
zone  width,  it  is  a  reasonable  approximation  to  treat  the  re¬ 
action  as  instantaneous.  However,  for  a  wide  variety  of  ex¬ 
plosives,  including  blasting  slurries,  ANFO,  and  other  composite 
explosives,  this  is  not  true.  For  a  typical  fast-reacting 
explosive  (e.g.,  HMX ,  LX-04)  the  reaction  time  is  on  the  order 
of  0.1  usee  whereas  for  a  typical  ANFO  mixture  (94  percent 
prilled  ammonium  nitrate  to  6  percent  No.  2  diesel  fuel  oil)  the 
reaction  time  is  on  the  order  of  100  sec  or  greater.  In  this 
case,  for  most  charge  sizes  of  interest,  reaction  rate  effects 
must  be  taken  into  account. 

3.2  REACTION  RATE  MODELS 

One  of  the  simplest  models  for  explosive  reactions  is  the 
Arrhenius  first-order  reaction  rate.  The  first-order  rate 
equation  is  (Reference  7) 


where  F  denotes  the  mass  fraction  of  the  constituents  that  have 
reacted,  and  hence  also  the  fraction  of  the  total  available  heat 
of  reaction  released.  The  Arrhenius  expression  for  the  rate 
constant  is 


-B/BT 
k  -  Ae  a 


(17) 


where  R  is  the  universal  gas  constant,  T  is  the  temperature , 
is  tne  activation  energy,  and  A  is  the  frequency  factor  for  the 
reaction.  The  factor  A  may  depend  weakly  on  temperature  accord¬ 
ing  to  the  relation 


_  __a 

A  »  BT 


(18) 


where  B  *  constant  and  0  <  a  H  1.  Here  we  will  consider  A 
independent  of  temperature  and  rewrite  (17)  as 


k 


1 

T 


•E  /RT 


(19) 


where  x  *  1/A  is  the  characteristic  reaction  time. 

The  activation  energy  represents  a  potential  energy  barrier 
for  the  reaction.  This  is  shown  schematically  in  Figure  18  from 
Reference  7.  When  the  thermal  energy  of  the  reactants  is  small 
compared  to  E  ,  the  exponential  term  in  the  rate  constant  is 
near  zero  and  the  reaction  will  not  proceed.  When  the  tempera¬ 
ture  in  the  explosive  increases,  due  to  shock  heating,  the  rate 
constant  increases  and  the  reaction  proceeds. 


The  first-order  rate  equation  (Equation  16)  is  primarily 
applicable  to  gaseous,  liquid,  or  homogeneous  solid  explosives, 
i.e.,  explosives  in  which  the  constituents  are  available  for 
reaction  instantaneously  following  the  initial  shock  heating. 
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For  ANFO,  however,  the  unreacted  explosive  is  in  the  form  of 
porous  ammonium  nitrate  prills,  approximately  1  to  2  mm  in 
diameter,  coated  with  fuel  oil.  For  the  constituents  to  react 
the  ammonium  nitrate  prills  must  first  decompose  and  then  mix 
with  the  fuel  oil.  These  processes  are  not  taken  into  account 
in  the  first-order  rate  equation. 

Several  models  have  been  proposed  to  describe  the  reaction 
process  in  an  ANFO  detonation.  If  one  assumes  that  the  decompo¬ 
sition  of  the  ammonium  nitrate  prill  is  the  controlling  process 
and  that  the  mixing  of  the  decomposition  products  and  the  fuel 
oil  takes  place  essentially  instantaneously  as  the  decomposition 
proceeds,  then  the  grain  burning  model  proposed  by  Eyring  (Ref¬ 
erence  8)  is  a  reasonable  model.  On  the  other  hand.  Finger 
(Reference  9)  believes  that  the  decomposition  takes  place  quickly 
and  that  the  reaction  rate  is  controlled  by  the  diffusion  of  the 
reactants.  In  this  work  we  have  considered  only  the  grain  burn¬ 
ing  model. 

In  the  grain  burning  theory,  the  rate  equation  is  dependent 
on  the  shape  of  the  grains  or  prills  and  on  the  type  of  ignition. 
For  spherical  grains  ignited  uniformly  over  the  surface,  the 
rate  equation  is  (Reference  10) 

=  3k (1  -  F)2/3  (20) 

where  k  is  given  again  by  Equation  19.  Under  isothermal  condi¬ 
tions  this  may  be  integrated  to  give 


where  x  is  the  time  for  complete  reaction.  Figure  19  from 
Reference  11  shows  examples  of  other  configura cions,  i.e.,  a 
sphere  ignited  at  the  center,  a  sphere  ignited  at  one  point  on 
the  surface,  and  a  plane  ignited  on  two  sides.  It  is  clear  that 
the  amount  of  reacted  material  at  any  given  time  increases  with 
the  number  of  ignition  points.  Since  these  ignition  points 
typically  occur  at  points  of  contact  between  grains,  for  a  mix¬ 
ture  of  spherical  prills  of  ANFO,  ignition  should  occur  at  several 
points  (typically  6  to  12  according  to  Reference  11)  on  the  sur¬ 
face.  For  this  case  a  uniform  ignition  over  the  surface  is  not 
a  bad  approximation. 

There  is  a  basic  difference  between  the  grain  burning 
theory  and  the  first-order  rate  equation.  Grains  burn  completely 
in  the  time  x.  According  to  the  first-order  rate  equation, 
however,  the  fraction  of  reacted  constituents  increases  exponen¬ 
tially  to  1  with  a  cx»sracteristic  time  x.  A  comparison  of  the 
two  models  is  shown  in  qure  20. 

To  evaluate  the  effect  of  the  activation  energy  and  the 
reaction  time  constant  for  both  this  first-order  rate  model  and 
the  spherical  grain  burning  model,  a  series  of  finite-difference 
calculations  was  performed  with  the  one-dimensional  POD  code. 
Finite-difference  representations  of  Equations  16  and  20  were 
programmed  into  the  POD  code  with  the  rate  constant,  k,  given  by 
an  equation  equivalent  to  Equation  19.  For  convenience,  the 
activation  energy  E  in  Equation  19  was  defined  by  a  dimension- 
less  fraction,  a,  according  to 

E  =  otQ  (22) 

cl 
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Figure  19  Burning  rate  for  various  models  (from  Reference  11). 
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Figure  20  Comparison  of  isothermal  burning  for  various 
models . 
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and  Equation  19  was  rewritten, 


k 


1  e-otQ/E 
T 


(23) 


where,  again  for  convenience  in  programming ,  the  factor  RT  has 
been  replaced  by  the  internal  energy  E. 


A  series  of  one-dimensional  (plane  symmetry)  calculations 
was  made  with  the  first-order  rate  model  to  evaluate  the  effects 
of  variations  in  the  activation  energy  parameter,  a,  and  the 
time  constant,  x.  In  all  cases,  the  steady-state  detonation 
properties  were  the  same;  so  the  *inal  values  of  detonation 
velocity  should  be  the  same.  Detonation  front  trajectories  for 
several  values  of  a  are  shown  in  Figure  21.  This  figure  shows 
that  the  va) ue  a  is  most  important  to  the  initiation  of  the 
detonation  and  has  little  effect  once  the  full  detonation  wave 
has  developed.  As  a  increases  the  detonation  builds  up  more 
slowly  until  for  a  ~  0.5  detonations  could  not  be  initiated. 
Values  of  a  between  0.01  and  0.1  gave  good  results  for  calcula¬ 
tions  of  ANFO  detonations. 


Figure  22  shows  detonation  trajectories  computed  with 
different  time  constants  in  the  first-order  rate  equation.  As 
expected,  long  reaction  times  result  in  a  slow  buildup  of  the 
detonation  wave.  For  x  =  100  psec,  steady  state  was  never 
achieved  within  the  propagation  distance  calculated. 

A  similar  series  of  calculations  was  performed  using  the 
grain-burning  model.  Figure  23  shows  the  effects  of  variations 
in  x.  The  general  effect  is  the  same  as  for  the  first-order  rate 
model.  By  adjusting  x,  it  was  possible  to  fit  an  experimentally 
determined  trajectory  for  undiluted  ANFO  from  one  of  the  small- 
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Figure  21  Evaluation  of  the  effect  of  activation  energy  on 
calculated  detonation  trajectories  for  the  first- 
order  reaction  rate  model  (x  held  constant  =  20  usee) . 
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Figure  22  Evaluation  of  the  effect  of  the  reaction  time  constant 
t,  on  calculated  detonation  trajectories  for  the  first 
order  reaction  rate  model  (a  held  constant  =  0.1). 
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Figure  23  Evaluation  of  the  effect  of  the  reaction  time  constant, 
t,  on  calculated  detonation  trajectories  for  the  grain 
burning  model  (a  held  constant  =  0.1). 
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scale  tests  using  this  spherical  grain-burning  model.  This  is 
shown  in  Figure  24.  More  comparisons  with  experiment  will  be 
discussed  later. 

Typical  calculated  results  of  detonation  wave  buildup  for 
the  two  rate  models  are  shown  in  Figure  25.  The  pressure  builds 
up  as  the  detonation  approaches  steady-state  until  it  reaches  a 
value  higher  them  the  Chapman-Jouguet  pressure  (PCJ) .  This  peak 
corresponds  to  the  Von  Neumann  spike  discussed  previously.  The 
fraction  of  material  reacted  is  indicated  by  the  dotted  line 
(related  to  the  scale  on  the  right) .  As  the  reaction  goes  to 
completion  the  pressure  drops  until  it  reaches  PCJ  at  the  time 
that  F  reaches  1.  The  calculated  spike  pressure  is  not  physically 
realistic  in  these  calculations.  Its  magnitude  is  dependent  on 
the  Hugoniot  of  the  unreacted  explosive  (see  Figure  17) ,  which 
was  not  accurately  modeled,  and  on  dissipative  effects  in  the 
shock  front  which  were  not  included  here.  The  Chapman-Jouguet 
pressure  is  probably  a  better  estimate  of  the  actual  peak  deto¬ 
nation  pressure. 
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Figure  24  Comparison  of  measured  detonation  trajectory  with 
calculation  using  grain  burning  model. 
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Figure  25  Calculated  detonation  wave  profiles  for  two  reaction 
rate  models. 
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SECTION  4 


DETONATION  PROPERTIES  OF 
AMMONIUM  NITRATE/FUEL  OIL  EXPLOSIVES 

Ammonium  nitrate  fuel  oil  mixtures  have  been  widely  used  as 
blasting  agents  for  many  years  in  the  mining  industry  and  con-* 
siderable  research  has  been  performed  to  determine  their  detona¬ 
tion  characteristics.  This  work  is  summarized  in  References  10 
through  14.  While  this  work  was  in  general  quite  well  done,  *-he 
investigators  were  hampered  by  not  being  able  to  experiment  with 
charge  sizes  large  enough  to  eliminate  the  effects  of  the  long 
reaction  time.  Consequently,  detonation  velocities  and  pressure 
have  been  predicted  that  are  substantially  lower  than  those 
measured  for  large  charges. 

Recently  a  number  of  tests  have  been  conducted  using  large 
charges  of  ANFO.  The  Naval  Ordnance  Laboratory  (References  15 
and  16)  conducted  tests  to  determine  the  airblast  characteristics 
of  ANFO  compared  to  TNT;  a  series  of  small  tests  and  three  large 
tests (20  to  100  tons)  were  conducted.  Measurements  of  the  aver¬ 
age  detonation  velocity  of  the  charges  were  made  which  agreed 
well  with  values  reported  in  the  literature  cited  above.  Un¬ 
fortunately,  no  detailed  measurements  of  the  detonation  trajectory 
(or  detonation  velocity  as  a  function  of  distance  of  propagation) 
were  made.  On  the  MINE  THROW  I  event  (120  tons)  detonation 
velocities  and  pressures  were  measured  which  were  much  higher 
than  previously  observed  (Reference  17) .  The  detonation  velocity 
continued  to  increase  over  the  entire  distance  of  propagation, 
reaching  a  final  value  of  ~  0.6  cm/ysec.  Pressures  measured  by 
manganin  gauges  on  the  outer  surface  of  the  charge  were  on  the 
order  of  90  to  100  kbar. 


PrmAv  pift  Milk 
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Although  these  date  disagree  with  previous  measurements, 
they  are  the  first  detailed  detonation  velocity  and  pressure 
measurements  made  on  large  ANFO  charges  and  probably  are  more 
representative  of  the  steady-state  behavior  of  ANFO.  This  is 
supported  by  data  from  the  test  series  described  in  this  report 
and  by  predictions  of  ANFO  detonation  properties  based  on  an 
empirical  method  developed  by  Kamlet  et  al.  (References  1 
through  4)  . 

4.1  EMPIRICAL  PREDICTIONS 

Kamlet  and  his  co-workers  have  developed  a  simple  empirical 

method  for  predicting  the  detonation  properties  of  C-H-N-0 

explosives.  They  have  demonstrated  that  for  a  wide  variety  of 

3 

common  military  explosives  with  pQ  >  1.0  gm/cm  the  method 
agrees  well  with  data  and  with  the  predictions  of  the  RUBY  code. 
Also,  Finger  (Reference  18)  has  found  that  the  method  agrees 
well  with  his  data  on  slurry  explosives.  The  empirical  equations 
are 


P 


(24) 


D  =  Ag> **  (1  +  B  pQ) 
<p  =  NM”  (T 


(25) 

(26) 


where  K  =  15.58,  A  =  1.01,  B  =  1.30,  P  is  the  detonation  pressure 
in  kilobars,  D  is  the  detonation  velocity  in  millimeters/psec, 

N  is  the  moles  of  gas  per  gram  of  explosive,  M  is  the  average 
molecular  weight  of  the  detonation  product  gases  in  grams  of  gas 
per  mole  of  gas,  Q  =  -AHQ  is  the  chemical  energy  of  the  detona¬ 
tion  reaction  in  calories  per  gram,  and  pQ  is  the  initial  density 
of  the  explosive  in  grams  per  cubic  centimeter. 
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Although  these  equations  are  generally  applicable  only  for 
densities  above  1.0  gm/cm3,  we  have  attempted  to  apply  them  to 
ANFO  at  lower  densities.  The  density  restriction  comes  from  an 
assumption  about  the  reactions 


2C0 

4- 

co2 

+  c, 

AH  *-41.2  kcal 
o 

(27) 

H2 

-► 

4- 

h20 

+  c, 

AH  *-31.4  kcal 
o 

(28) 

These  are  assumed  to  be  predominantly  to  the  right  which  is  true 
only  at  higher  pressure  (higher  loading  densities).  However,  the 
reaction  considered  here  (Equation  29  below)  is  oxygen  balanced 
so  that  no  free  carbon  exists  in  the  products  and  the  above 
reactions  do  not  occur.  Therefore,  the  equations  should  be  valid 
for  ANFO  at  lower  densities. 

To  calculate  the  detonation  properties  of  ANFO  we  assume 
the  following  reaction  (for  a  stoichemetric  mixture  of  94  per¬ 
cent  ammonium  nitrate  and  6  percent  No.  2  diesel  fuel  oil) : 

3NH4N03  +  CH2  ♦  3N2  +  7H20  +  C02  (29) 

For  this  reaction  N  *  1/M  =  0.0433  mcles/gm,  Q  =  912  cal/gm. 

The  resultant  equations  for  detonation  velocity  and  pressure  are 

D(mm/Msec)  *  3.292  pQ  +  2.532  (30) 

P  =  97.95  pQ2  (31) 

which  are  plotted  in  Figure  26  for  the  range  of  densities  of 
interest  here. 
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The  Chapman- Jouguet  pressure  (calculated  by  solving 
Equations  10  and  11  for  P  as  a  function  of  D,  assuming  Q  *  912.8 
c&i/gm)  corresponding  to  the  detonation  velocity  determined  by 
Equation  30  is  shown  as  a  dashed  line  for  comparison.  The 
empirical  equation  predicts  slightly  lower  detonation  pressures 
than  would  be  predicted  by  Chapman- Jouguet  theory. 

4.2  COMPARISON  WITH  EXPERIMENT 

Figure  26  also  shows  a  comparison  of  the  predictions  with 
data  from  the  series  of  large-scale  tests  of  diluted  ANFO.  The 
detonation  velocities  are  the  final  values  after  steady-state  had 
been  achieved,  and  the  pressures  are  the  calculated  Chapman- 
Jouguet  pressures  corresponding  to  those  detonation  velocities. 
The  measured  detonation  velocities  are  on  the  order  of  10  percent 
lower  than  the  predicted  values  while  the  pressures  are  generally 
in  slightly  better  agreement. 

The  agreement  between  measured  and  calculated  detonation 
velocities  is  not  as  good  as  was  reported  in  Reference  4  for 
other  explosives.  There  are  several  factors  that  could  be 
responsible  for  this.  First  of  all,  our  contention  that  the 
technique  can  be  used  for  pQ  <  1.0  gm/cnr  for  ANFO  may  be  in¬ 
correct;  perhaps  there  are  other  factors  not  considered  that 
make  the  technique  invalid  in  this  range.  Secondly,  the  ex¬ 
plosives  to  which  the  technique  was  applied  in  References  1 
through  4  were  generally  fine-grained,  homogeneous  explosives 
that  are  apt  to  behave  somewhat  differently  from  a  more  hetero¬ 
geneous  explosive  such  as  ANFO.  Finally,  it  is  possible  that 
the  data  obtained  from  the  large-scale  tests  did  not  represent 
steady-state  values  and  that  for  larger  charge  sizes  (longer 
distances  of  propagation)  higher  detonation  velocities  would  be 
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found.  The  irregularity  in  the  measured  detonation  trajectories 
(presumably  caused  by  incomplete  mixing  within  the  bags)  makes  it 
impossiblo  to  be  absolutely  certain  that  a  constant  velocity  has 
been  achieved.  However,  in  developing  a  model  for  the  detonation 
properties  we  will  assume  that  steady-state  has  been  reached  and 
that  the  data  from  the  tests  is  a  better  representation  of  the 
steady-state  behavior  of  ANFO  than  the  empirical  predictions. 

4.3  A  DETONATION  MODEL  FOR  ANFO 

From  the  test  data  a  mathematical  model  for  the  detonation 

properties  of  diluted  ANFO  was  developed.  A  density  of  0.75 

3 

gm/cm  was  selected  for  use  in  the  MINE  THROW  II  Event.  Data  from 
the  large-scale  tests  indicate  that  this  should  have  a  detonation 
velocity  of  0.47  cm/ysec.  A  polytropic  gas  (constant  y)  descrip¬ 
tion  of  the  behavior  of  the  detonation  products  was  used  since 
no  data  exist  to  justify  a  more  complicated  model.  The  variables 
at  the  Chapman- Jouguet  state  for  this  explosive  (calculated  from 
Equations  10  through  14)  are  given  in  Table  1. 

A  series  of  one-dimensional  (planar)  finite-difference 
calculations  was  performed  to  determine  which  rate  model  gave 
the  best  fit  to  the  measured  detonation  trajectories.  The  grain¬ 
burning  model  with  the  parameters  given  in  Table  1  gave  best 
agreement.  A  comparison  between  calculated  and  measured  detona¬ 
tion  trajectories  is  shown  in  Figure  27.  The  agreement  is 
reasonably  good,  although  with  small  changes  in  the  model 
parameters  a  better  fit  could  probably  be  achieved. 
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Time ,  usee 


Distance,  cm 


Figure  27  Comparison  of  measured  and  calculated  detonation 
trajectories  for  diluted  ANFO. 
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TABLE  1 


MODEL  PARAMETERS  FOR  DILUTED  ANFO 


3 

pQ  =  0.75  gm/cm 

D  =  0.47  cm/ysec 

Q  =  912.8  cal/gm  =  0.03821  Mbar-cm'Vgm 


Chapman- Jouguet 
State  Variables 


Reaction  Rate  Model 
Parameters  (Grain-Burning  Model) 


CJ 

YCJ 

UCJ 

vCJ/vo 


=  0.0558  Mbar 
■  1.97 

=  0.158  cm/ysec 
*  0.312 


a  =  0.05 
x  =  100  ysec 


Although  this  model  gives  good  agreement  with  the  data,  this 
does  not  mean  that  it  necessarily  represents  a  physically  correct 
model  of  the  ANFO  detonation  process.  Other  models  for  the 
reaction  process  and  for  the  equation  of  state  of  the  reaction 
products  may  be  more  physically  realistic.  However,  the  objective 
of  this  program  was  to  find  a  simple  model  that  could  adequately 
describe  the  ANFO  detonation  properties  in  finite  difference 
calculations  of  the  MINE  THROW  II  simulation  charge  design. 

Further  modeling  work  is  beyond  the  scope  of  this  program  but  is 
recommended  for  future  study. 
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SECTION  5 

SUMMARY  AND  CONCLUSIONS 


The  test  series  discussed  in  Section  2  demonstrated  the 
feasibility  of  producing  a  low-detonation-velocity  ANFO  explo¬ 
sive  by  diluting  normal-density  ANFO  with  a  low-density  inert 
material.  Polystyrene  beads  were  found  to  be  the  best  diluent 
from  the  point  of  view  of  minimizing  both  cost  and  mass  of 
inert  material  added  to  the  explosive.  Although  the  beads  tend 
to  segregate  from  the  ANFO,  by  mixing  the  correct  proportions 
of  beads  and  ANFO  in  each  bag,  the  scale  of  the  density  varia¬ 
tions  can  be  minimized.  Some  irregularities  in  the  detonation 
trajectories  of  the  low-density  mixtures  were  observed  which 
were  apparently  caused  by  incomplete  mixing  within  the  bags. 

To  prevent  possible  problems,  a  better  mixing  procedure  should 
be  developed  for  the  explosive  produced  for  MINE  THROW  II. 

A  mixture  having  a  density  of  0.75  gm/cm  was  selected  for 
MINE  THROW  II.  The  detonation  pressure  of  this  explosive  is 
predicted  to  be  on  the  order  of  56  kilobars,  which  satisfies  the 
design  requirements  for  the  simulation  charge.  Furthermore,  less 
than  1  percent  by  mass  of  the  polystyrene  diluent  is  required  to 
achieve  this  density;  this  minimizes  any  effect  the  diluent  may 
have  on  the  detonation  reaction  or  on  the  behavior  of  the 
reaction  products. 
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A  mathematical  model  of  the  detonation  properties  of  the 
low-density  ANFO  was  developed.  A  Chapman- Jouguet  detonation 
was  assumed  and  a  grain-burning  reaction  rate  model  (Reference  8) 
war  used  to  model  the  reaction  process.  The  data  indicated 
that  a  distance  of  propagation  of  100  to  150  cm  was  required  to 
establish  a  steady-state  detonation  which,  for  the  reaction  rate 
model  used,  corresponds  to  a  reaction  time  of  100  usee.  The 
detonation  parameters  of  the  explosive  to  be  used  for  the 
MINE  THROW  II  Event  are  presented  in  Table  1. 
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